Dispersive and Covalent Interactions Between Graphene and Metal Surfaces from the 

Random Phase Approximation 
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We calculate the potential energy surfaces for graphene adsorbed on Cu(lll), Ni(lll), and 
Co(0001) using density functional theory and the Random Phase Approximation (RPA). For these 
adsorption systems covalent and dispersive interactions are equally important and while commonly 
used approximations for exchange-correlation functionals give inadequate descriptions of either van 
der Waals or chemical bonds, RPA accounts accurately for both. It is found that the adsorp- 
tion is a delicate competition between a weak chemisorption minimum close to the surface and a 
physisorption minimum further from the surface. 

PACS numbers: 71.15.Nc, 73.22.Pr, 81.05.ue 



The recent experimental realization and characteriza- 
tion of isolated graphene sheets [2, Q , have boosted a ma- 
jor interest in this novel two-dimensional material. The 
most prominent feature is perhaps the band structure, 
which exhibits linear dispersion near the Fermi level re- 
sulting in a relativistic description of the charge-carrying 
quasi-particles In particular, graphene shows a re- 
markably high intrinsic carrier mobility, and therefore 
seems very well suited for nanoscale electronics devices. 
For such applications, the coupling to metal contacts 
plays a fundamental role and measurements show that 
graphene binds very differently on various metal sur- 
faces. Understanding the interactions between graphene 
and metal surfaces, therefore becomes a most important 
task since the adsorption geometry and bond distance 
may have drastic consequences for the electronic struc- 
ture and transport properties of isolated graphene layers. 
For example, Pd(lll), Co(0001), and Ni(lll) have been 
demonstrated to induce a band gap in adsorbed graphene 
sheets, which signals a covalent bond with the metal 
In contrast, adsorption on Cu(lll), Ag(lll), 
Au(lll), and Pt(lll) do not change the graphene band 
structure significantly [7H9(. Furthermore, graphene rep- 
resents a prototypical example of a 7r-conjugated sys- 
tem. Unraveling the coupling mechanism to different 
metal surfaces, could therefore possibly improve the un- 
derstanding of adsorption mechanisms for a large variety 
of organic molecules on metals. 

Previous attempts to simulate the interaction between 
metal surfaces and graphene have employed Density 
Functional Theory (DFT) with various app roximations 
for the exchange-correlation functional [id [Til ] . While 
the Local Density Approximation (LDA), reproduces the 
strong binding at Pd(lll), Co(0001), and Ni(lll) sur- 
faces it also predicts a strong coupling to Cu(lll), which 
has been shown to bind graphene weakly [8j. Further- 
more, the Generalized Gradient Approximations (GGA), 
which supposedly improves the LDA description, gives 
significantly different results than the local density ap- 



proximation, and the LDA bonding therefore seems to 
be fortuitous rather than contain the essential physics. 
In this respect, the major shortcoming of the semi-local 
functionals is the lack of van der Waals interactions, 
which are expected to be important for the adsorption 
of graphene on metals. On the other hand, applying the 
celebrated van der Waals functional (vdW-DF) [12| re- 
sults in a shallow physisorption minimum at all metal 
surfaces 1(1 U|, which can be traced to the inaccurate 
description of covalent bonds with this functional. Ad- 
sorption of graphene at metal surfaces thus represents 
an electronic structure problem where both LDA, GGAs, 
and the vdW-DF fail due to the detailed balance between 
covalent and dispersive interactions. 

In this letter we go beyond the semi-local and vdW-DF 
approximations for exchange and correlation and com- 
bine exact exchange (EXX) with a correlation energy ob- 
tained from the Random Phase Approximation (RPA). 
This functional has been shown to give an accurate de- 
scription of both van der Waals interactions and covalent 
bonds Il3l4la| as well as adsorption of molecules on metal 
surfaces [1 7| . We obtain the binding energy and bond dis- 
tance of graphene adsorbed on Cu(lll), Co(0001), and 
Ni(lll) and show that the results are in good agreement 
with experiments. 

Using the adiabatic connection and fluctuation- 
dissipation theorem (ACDF), the correlation energy can 
be evaluated in the Random Phase Approximation as: 
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X KS = xf S +xf S is the Kohn-Sham response 
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function and v is the coulomb interaction. To get total 
energies we combine the RPA correlation energy with the 
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exact exchange energy, which we obtain from 
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where and f/wkfr) are Kohn-Sham eigenvalues and 
eigenfunctions respectively and f{e ari k) are the occupa- 
tion numbers for ?/Wk(i")- Eq. ([2]) is derived from the 
ACFD and differs from the standard expression for ex- 
act exchange energy, if the occupation numbers are not 
integer valued. However, as discussed in Ref. [ltjj ]. it is 
natural to apply Eq. ^ when the exact exchange energy 
is combined with the RPA correlation energy, since the 
two expressions go hand in hand through the derivation 
of Eq. (JT]) . It has also been shown empirically that Eq. 
([2]) is less sensitive to the width of the artificial smearing 
function /(e) than the standard expression for exact ex- 
change. Furthermore, for metals the integrands in both 
expressions above diverge in the limit of q — > 0, and it 
has been demonstrated that if the q = terms in both 
expressions are excluded, the sum of Eqs. ((TJ)- ((2J) ex- 
hibits fast convergence with respect to k-point sampling 
[l6j . The frequency integration in Eq. ([1} is carried out 
using 16 Gauss-Legendre points with a weight function 
ensuring that the integral of f(x) oc x ( 1 / B - 1 ) cx P~ axl/B 
is reproduced exactly. We have used B = 2.5 and a 
is determined by the frequency cutoff, which we set to 
800 eV. With this frequency sampling the RPA correla- 
tion energy is converged to within a few meV. Since, the 
present approach is not self-consistent, one has to choose 
a set of orbitals, on which Eqs. CO)-© are evaluated and 
for all the calculations below we have used self-consistent 
PBE orbitals. We have compared the RPA potential en- 
ergy surface for graphene on Ni(lll) using self-consistent 
PBE orbitals with the result obtained with self-consistent 
LDA orbitals and the results are nearly indistinguishable. 
In the following we will refer to the Hartree-Fock energy 
E HF = Epbe - E£ C BE + E® xx , where E PBE is the re- 
sult of a self-consistent PBE calculation and E EXX is 
evaluated on the self-consistent PBE orbitals. 



The calculations were performed with gpaw 1 
which is a Density Functional Theory code using the 
projector augmented wave method [2~l| and uniform real- 
space grids. The response function and Hartree Fock en- 
ergy were calculated in a plane wave basis after having 



Fourier transformed the real space wavefunctions 22 1. 
All self-consistent calculations were carried out with a 
grid spacing of 0.18 A. For the exact exchange energy 
we used a plane wave cutoff of 870 eV and for the re- 
sponse function we used a cutoff of 150 eV. The number 
of bands included in the response function were set equal 
to the number of plane waves defined by the cutoff en- 
ergy. We tested a few RPA energy differences using a 
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FIG. 1: (color online). Potential energy surfaces of graphene 
on Cu(lll), Ni(lll), and Co(0001). The adsorption geometry 
for graphene on fcc(lll) is shown in the top left. 



200 eV cutoff and found that the results were converged 
to within ~ 2 meV per carbon atom. The metal surface 
was simulated using four layers where the two top layers 
were relaxed using the PBE functional and we checked 
that the results obtained with LDA and vdW-DF do not 
change if the slab thickness is increased to six layers. The 
graphene lattice constant (a = 2.46 A) were scaled to fit 
the experimental minimal surface unit cell of the met- 
als, which have side length of 2.49 A, 2.51 A and 2.56 
A for Ni, Co, and Cu respectively. Due to the nonlocal 
nature of the correlation energy a large amount of vac- 
uum is needed in the non-periodic direction and we found 
that the calculations are converged when the metal slabs 
were separated by 20 A of vacuum. For the Ni(lll) and 
Co(0001) slabs the calculations were spin-polarized. A 
12 x 12 gamma-centered k-point mesh was used for Ni 
and Co whereas an 8 x 8 grid was sufficient for Cu, but 
we return to the issue of k-point convergence below. 

The potential energy surfaces for graphene on Cu(lll), 
Co(0001), and Ni(lll) are shown in Fig. HJ The C atoms 
were fixed at a top and a fee hollow site for Cu and Ni 
and top and hep hollow for Co. vdW-DF gives almost 
identical binding for the three metals with a minimum 
at d = 3.7 A and a binding energy of ~ 40 meV per C 
atom. LDA shows strong binding at d ~ 2.0 A on all 
three surfaces and PBE gives rise to weak binding at the 
Co(0001) surface only. For the Cu(lll) surface, RPA 
predicts a physisorption minimum at d ~ 3.25 A with 
a binding energy of 62 meV per C atom. While there 
is no direct experimental measurements for graphene on 
pure Cu(lll), the result is in good agreement with ex- 
perimental studies where Cu was intercalated between 
graphene adsorbed on Ni(lll), resulting in the graphene 
sheet becoming weakly coupled to the metal [8J, [23 1 . In 
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FIG. 2: (color online), fc-point dependence of the energy 
difference between the chemisorbed states (d = 2.25 A) and 
the physisorbed state (d = 3.0 A) for graphene on Ni(lll). 



contrast, graphene has been found to adsorb on both 
Ni(lll) and Co(0001) in registry with the minimal sur- 
face unit cell. This is due to the close match of lattice pa- 
rameters between graphene and these surfaces, and PBE 
calculations confirm that the stretching energy needed 
for graphene to match the minimal unit cell is smaller 
than the calculated binding energy. In Co and Ni the 
d-band straddles the Fermi level and hybridization with 
the graphene 7r-bands makes it possible to form chemi- 
cal bonds at these surfaces. This is indeed what we find 
with the RPA functional. A physisorption minimum at 
d = 3.25 A is still observed at both surfaces, but there is 
also a competing chemical interaction closer to the sur- 
face. The RPA functional captures both effects and the 
calculated potential energy surface can be regarded as a 
result of these competing binding mechanisms. At the 
Co(000f) surface, RPA predicts a global minimum at 
d = 2.3 A, which is in good agreement with the experi- 
mental binding distance The minimum represents a 
chemical bound state with a binding energy of 86 meV 
per C atom. The physisorbed state at d = 3.25 A is ob- 
served as a weak local minimum at d = 3.25 A and is 
situated 12 meV higher than the chemisorbed state. In 
the case of Ni(lll) the physisorbed state at d = 3.25 A 
has a slightly lower energy 8 meV than the chemisorbed 
state at d = 2.3 A. The experimental binding distance 
for graphene on Ni(lll) ranges from d = 2.1 A [24| to 
d = 2.8 A [25| and while band structure measurements 
indicate that the value is probably closer to the former 
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111 ], the disagreement could suggest a broad flat 
minimum in the potential energy surface as predicted by 
the RPA functional. Finally, it should be noted that the 
RPA potential energy surface matches the vdW result for 
large distances in all three cases, which demonstrates the 
excellence of vdW-DF when dealing with pure van der 
Waals forces. 

Since the binding energy of the physisorbed and 
chemisorbed states at Ni(lll) are only marginally sepa- 
rated in energy, we have checked if a denser fc-point sam- 
pling could change the preferred binding distance. The 
energy difference E(d = 2.25 A) - E(d = 3.0 A) is shown 
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FIG. 3: (color online). Left: RPA correlation contributions 
to the total RPA energy. Right: HF contribution to the total 
RPA energy. 



as a function of fc-point sampling in Fig. [21 The result is 
seen to be well converged at 16 x 16 fc-points where the 
energy difference is 4 meV (in favor of the physisorbed 
state). However, such a small energy difference is be- 
yond the accuracy of the RPA functional and it cannot 
be concluded that the physisorbed state is the preferred 
one. Fig. [5] also nicely confirms the rapid convergence 
of Eh f + E^ PA , despite the fact that either functional 
exhibits rather slow convergence. 

The expected chemisorption minimum at d = 2.1 A 
for Ni(lll) |24j is not reproduced by RPA. However, it is 
well known that RPA tends to underestimate atomization 
energies for molecules [26} and this trend has also been 
observed for adsorption of CO on various metal surfaces 
17| . On the other hand, RPA accurately describes the 



van der Waals interaction between graphene layers (14 1. 
One could thus expect that the physisorption minimum 
for graphene on metals is rather well described by RPA, 
whereas the chemisorption minimum is most likely under- 
estimated. This deficiency has been shown to be partly 
cured by inclusion of second order screened exchange [27j | 
and self-consistency corrections [28| . In the present case, 
a correction of the RPA underbinding would therefore 
most likely lead to the energy of the chemisorption min- 
imum being lowered compared to the physisorbed state. 
For Ni(lll), this would change the global minimum to 
the chemisorbed state. 

Contrary to the physisorbed state at Cu(lll), the 
chemical interaction for graphene on Ni(lll) and 
Co(000T) are expected to be highly site dependent and 
this is indeed observed. The RPA potential energy sur- 
face has been calculated for Ni(lll), where the C atoms 
were put at the fee and hep hollow sites instead of a top 
and fee hollow site. The result then resembles the RPA 
potential energy surface for graphene on Cu(lll) with a 
distinct physisorption minimum at d — 3.25 A indicating 
that the top carbon atom gives rise to the covalent bond. 

Although the RPA potential energy surfaces for 
graphene on the three surfaces seem rather similar, 
there arc major qualitative differences in the mechanism 
that produces binding, and the contributions from the 
Hartree-Fock energy and the RPA correlation are very 
different in the three cases. In Fig. [3] we show the sep- 
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FIG. 4: (color online). Band structure of graphene the three 
metal surfaces. The dots represent the projected density of 
state of the ir and ir* bands of graphene. For Co and Ni the 
majority spin channel is shown. The minority spin channels 
have a similar structure. 



arate contributions from HF and RPA. In particular, it 
should be noted that while the potential energy surfaces 
for graphene on Ni(lll) and Co(0001) are very similar, 
the HF calculation produces a chemisorption minimum 
for Ni(lll) whereas HF is purely repulsive at Co(0001). 
Since the van der Waals interaction can be regarded as an 
attraction between density fluctuations, one could spec- 
ulate the differences are due to matching of the surface 
plasmon frequencies of the metals with that of graphene 

2.25 A have tenta- 



Until now, binding distances of d ■ 
tively been referred to as chemisorbed, while distances 
of d ~ 3.25 A have been referred to as physisorbed. 
We will now justify this interpretation in the sense that 
chemisorbed states show hybridization between graphene 
and metal bands while physisorbed states do not. In 
Fig. HI the projected density of states for the graphene 
7r and 7T* bands on the three metal surfaces are shown. 
It is observed that the graphene gap is opened at the 
chemisorbed state, whereas the band structure is nearly 
conserved at the physisorbed state. 

In conclusion, we have calculated the potential energy 
surfaces for graphene on Cu(lll), Ni(lll), and Co(0001) 
using the RPA method. The potential energy surfaces 
are very different from those obtained with local and van 
der Waals approximations for the exchange-correlation 
energy and clearly show the detailed balance between 
dispersive and chemical interactions. 
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